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Abstract -In one-channel, finite-size Luttinger one-dimensional quantum dots, both Friedel os¬ 
cillations and Wigner correlations induce oscillations in the electron density with the same wave¬ 
length, pinned at the same position. Therefore, observing such a property does not provide any 
hint about the formation of a Wigner molecule when electrons interact strongly and other tools 
must be employed to assess the formation of such correlated states. We compare here the behavior 
of three different correlation functions and demonstrate that the integrated two point correlation 
function, which represents the probability density of finding two particles at a given distance, is 
the only faithful estimator for the formation of a correlated Wigner molecule. 


Introduction. — Strongly interacting one¬ 
dimensional (ID) systems are attracting an increasing 
attention: ID quantum dots in carbon nanotubes [1,2 
semiconducting heterostructures and nano-wires 3, 4 
helical ID edge states in two dimensional topological 
insulators [5j-[8| and in spin-orbit coupled quantum 
wires |9 10 and the cold atom simulators li , are only 
some of the most prominent examples. 

From the experimental point of view, ID quantum 
dots (12] represent an invaluable tool: their small di¬ 
mension reinforce interaction effects, allowing for their 
detection. Moreover with transport spectroscopy one 
can investigate the low energy properties 


as the electron and spin density 13 -19 


From the theoretical perspective, ID interacting electrons 
have been intensely studied: aside from numerical tech¬ 
niques, such as exact diagonalization 20 21 , quantum 
Montecarlo 22 23, density functional theory 24-27 


and density matrix renormalization group 28 , exact 


solutions are possible thanks to the Bethe ansatz so¬ 
lutions 29 30 . Moreover a powerful held theory, the 

is also available. 


Luttinger liquid theory 31 33 


validity ranges from semiconducting quantum wires 34 
carbon nanotubes 35-38, edges in the integer and 


wires 


fractional quantum Hall effect 39j-j41 
topological insulators 

[9] [To 


to two dimensional 
6|[7], spin-orbit coupled quantum 
cold atoms in ID optical lattices 42 43 


spin systems 44 , and mesoscopic circuits 45 . The main 


feature of the Luttinger liquid is the bosonic character 
of its low energy excitations 31 , even for strong interac- 


12 , and, in 


combination with local probes, obtain information such 


tions. This issue allows for analytical calculations in any 
interaction regime. 

In order to describe finite-size, sharply confined systems, 
such as quantum dots, one has to to adopt finite-size 
boundary conditions (FBC), such as open-boundary 
conditions, which strongly differ from periodic boundary 
conditions (PBC) [32p7p6] , In general, the effect of FBC 
is to halve the number of channels of a Luttinger liquid 
theory with respect to the case of PBC. For example, 
finite-size spinful ID Luttinger liquids - which for PBC 
have four independent channels, exhibit only the two 
charge and spin channels, as reflections at the boundaries 
mixes left- and right-movers 32 . This system will be 


Its 


dubbed henceforth a two-channel Luttinger liquids (2LL). 
Similarly, systems sporting two independent channels with 
PBC will be enforced, by FBC, to become one-channel 
Luttinger liquids (ILL). A typical and very relevant 
example of this are the helical LL (HLL) occurring at the 
edge of topological insulators (5 -8 '. Here, spin-momentum 
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locking pairs the up (down) spin direction to right (left) 
movers yielding a 2LL with PBC. The additional presence 
of ferromagnetic barriers [47 ,[48] enforces twisted 
boundary conditions which mix spin and chirality of 
electrons 17 47 , resulting in a ILL. Other notable 


examples of ILLs are fully spin-polarized 2LLs and the 
spin-incoherent Luttinger liquid (SILL) 49-52 . The 


former can be easily achieved by applying a magnetic held 
to a 2LL. The SILL describes strongly interacting spinful 
electrons in the temperature regime D a <C ksT <C D p , 
D a and D p being the bandwidth of spin and charge 


excitations 53 . The corresponding Hamiltonian maps 
onto that of a ILL describing spinless fermions called 


holons 51 52 


Therefore, the ILL is far from being a pure theoretical 
model but represents nowadays a very interesting and 
lively subject of investigation. 

The most striking effect of strong, long range interactions 
in ID quantum dots is the formation of the Wigner 


molecule 54-58 , which is the finite-size counterpart of a 


Wigner crystal. Indeed, due to fluctuations, no Wigner 
crystal can occur in ID. However, when the electronic 
correlation length exceeds the length of the sample a 
correlated molecular state arises 55||56 . In a semiclassical 
picture, one can picture such a state as a regular array 
of N electrons, each of which free to oscillate around its 


equilibrium position 32 


In a 2LL, the density profile shows a competition between 
Friedel oscillations, a finite size effect present even in 
the absence of interactions, and Wigner oscillations, 
a pure interaction effect 15 24 28 58 [59 . The two 


oscillations are easily distinguishable on the basis of 
their different wavevectors, with Friedel (Wigner) oscil¬ 
lations being characterized by the wavevector fcp = 2fco 
(kw = 4fco)> with k 0 the Fermi momentum 151. It has 
been demonstrated that Friedel oscillations dominate for 
weak interactions, while Wigner oscillations dominate for 

However, while the presence 


32 33 


strong interactions 
of Wigner oscillations is a clear signature of interaction 
effects, the degree of correlation cannot be inspected 
within the electron density alone 27 . It has in fact been 


shown that in the presence of zero range interaction, 
Wigner oscillations appear in the electron density, while 
correlation functions unveil the uncorrelated nature of 


the state 27 


In ILL quantum dots on the other side both Friedel 
and Wigner oscillations have the same wavevector 
(2 ho) [49,|50}[60]. While in the noninteracting case it is 
easy to show, using the wavefunctions for a particle in 
a box, that the number of peaks of the local electron 
density of the system containing N particles is precisely 
N, in the case of strong interactions the N electrons of 
the Wigner molecule, oscillating around their equilibrium 
position, give rise again to N distinct peaks in the density. 
As a consequence it is not even possible to discriminate 
between finite size and interaction effects by studying 


the electron density alone 60 61 . The interplay between 


finite size and interaction effects and correlations must 
therefore be clarified in more details. 


The aim of this letter is to detect the transition from an in¬ 
teracting, liquid-like state to a strongly correlated Wigner 
molecule in a ILL. We identify a suitable tool to confirm 
the presence of a Wigner molecule in a ILL, comparing 
three different correlation functions and demonstrate that 
the probability density of finding two particles at a given 
distance is the best tool to detect Wigner correlations. 
By means of this, we show that strong interactions induce 
a crossover from an interacting but still liquid-like state 
to a Wigner molecule, in analogy to the case of a 2LL. We 
estimate the threshold value of the interaction parame¬ 
ter for the appearance of the Wigner molecule as g ~ 0.45. 

The outline of the Letter is the following: 

In order to elucidate the differences between 2LL and ILL, 
we start briefly reviewing the properties of the former, re¬ 
calling the crossover between Friedel and Wigner oscilla¬ 
tions in the electron density. Subsequently, we turn to the 
case of a ILL where we show that the electron density dis¬ 
plays no peculiar crossover when the interaction strength 
is increased. We then introduce three possible tools to de¬ 
tect the emergence of correlations and show that studying 
the conditioned probability to find two electrons at a given 
distance is among the most sensitive and reliable ways to 
assess Wigner correlations. 

Two-channel Luttinger liquid. — We start by re¬ 
calling the main results for a 2LL. The Hamiltonian H 2 LL 
reads [Ti = 1) [46] 



H 2 LL — H n + H b , 

(i) 

h n 

= Een 2 + ^N* 

2 2 

(2) 


OO 


H b 


(3) 


n—1 


Here, N = N + + N_, N a = N + — N_, N s (with s = ±) 
is the number of electrons with spin projection up/down 
and Hn represents the contribution of the zero modes, 
with E v = , Kv l ,/2Lg v written in terms of the velocity v v 
of the mode v = p, < 7 , of the system length L and of the 
Luttinger parameters g v . For repulsive interactions one 
has g p = g < 1, while g = 1 corresponds to the nonin¬ 
teracting limit. On the other hand, g a = 1 for an SU(2) 
invariant theory. The velocity v p = Vo/g of the charged 
mode (with Vo the Fermi velocity) is renormalized by the 
interactions while the spin modes velocity is v a = Vq. 

The term H b describes collective, quantized charge and 
spin density waves with boson operators d u ^ n and energy 
e v = 7 tv v /L. 
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The electron field operator satisfying open bound¬ 
ary conditions 'I , s (0) = = 0 is 46 4/ s (:r) = 

where ip St +(x) is a 2L-periodic fermion 


field that admits the bosonic representation 46 
V: 


ips,+( x ) = 


y/2- 


:e“ i9 - e i 


7ra 


• ®p(x) + s3>cr(x) 

e vi . (4) 


Here, a is the cutoff length, set as a = L/(ttN), 9 S satis¬ 
fies [0 S ,1V S '] = i6g iS i, and g s fulfill 7] s r) s ' + rjs'Vs = 2<5 SjS /, 
allowing the right anticommutation relations for different 
spins. The boson fields <J> p (cc), $ a (x) are given by 


OO _ mrac 


n—1 


nnx 

cos —-- ig v sm 


nnx \ 7+ 

— J <n+ h - c - 


(5) 

The particle density operator is p 2 ll(x) = ^ s=± p s (x) 
withp s (x) = 4 , J(a;)4' s (a;). Following a standard procedure 
it can be bosonized [3lj 58} 59,61,62; 


P2LL(x) = y + —d x <p p (x)+F Y P F S (x)+(l-F)p w (x). (6) 


Here 


Pf(x) 

p w (x) 

Vs(x) 
Pp/t t(x) 

f( x ) 



S=± 

N s 

2N s nx 

1 j 

L 

N \ 2Ntt 


— — COS 

1 j 


L 


- 2 f(x) - 2 ip s (x) 


- 4 f{x) - 2a/2 ip p (x) 


iPp(x) -(- Slfc,(x) 

V2 

^ [ < ^ > p/ct( — X ) ~ ^ > p/cr(2')] j 

1 i / sin(27ra:/L) 

- tan -1 -n—-- , , s 

2 y e ira/ L _ cos(2itx/L) 


,(7) 

,( 8 ) 

(9) 

( 10 ) 

( 11 ) 


where 0 < F < 1 models the relative weight of the Friedel 
and the Wigner contributions to the density [28] . 

From the above equations one can see that Friedel os¬ 
cillations display N/2 peaks for even N and N/2 + 1/2 
for odd N), while Wigner oscillations always display N 
peaks j59|. Figure [Tj shows the zero-temperature averaged 
electron density p 2L l{x) = (N, N<t\P 2 ll(x)\N, N„) for dif¬ 
ferent interactions. Here | N, N a ) is the ground state for N 
electrons and N a = 0 (N a = ±1) for even (odd) N. The 
average can be computed analytically 15 : results are not 
quoted here for simplicity. While in the non interacting 
case (g = 1) only Friedel oscillations appear, Fig. [lja), 
increasing interactions a competition between Friedel and 
Wigner oscillations develops as shown in Fig. 0b). Fi¬ 
nally, in the strong interaction regime (g —> 0) well devel¬ 
oped Wigner oscillations are present, see Fig. 0=)- Indeed, 
it can be analytically proven (not shown) that Wigner os¬ 
cillations dominate over Friedel ones in the density for 
9 < 1/3 [611. 

We can thus conclude that, even though the presence of a 
strongly correlated Wigner molecule must be inferred by 




Fig. 1: Plot of P2ll{x) (units 1 /L) as a function of x (units 
L) for N = 10 and different interaction strengths: (a) g — 1; 
(b) g = 0.4; (c) g = 0.1. In all panels F = 0.5, a = L/(10 tt), 
T= 0 . 


higher order correlation functions 27 , the crossover be¬ 


tween the weak and the strong interaction regimes can be 
inspected in terms of oscillations of the electron density. 

One-channel Luttinger liquids. — Let us now con¬ 
sider a ILL. For definiteness, we will consider the case 
of a fully spin-polarized 2LL although all the results pre¬ 
sented here apply to the other cases of ILL described in 
the introduction. The Hamiltonian is 


31 33 


Hill = + e 0 Y nb l b n, 


+oo 


( 12 ) 


n—1 


with Eq = nvo/Lg 2 and £ 0 = ttvq/L g. The electron 
operator 'F(x), satisfying 41(0) = 4/(L) = 0, is [46}|60] 
\F(a;) = ipR(x) — ipii(—x), with the 2L-periodic field ^r{x) 
given by 


i>R(x) = 


y/2i 


e -i6 e i^ e i<^(x) 


(13) 


Here the bosonic field <F(:r) is 


®(x) = Y 


n> 0 


e l 

Vgn L' 


/ nnx ' 


cos 


t 




L J 


ig sm (- 


/nnx\ 


L J\ 


b n + h.c. , 


(14) 

with [0,1V] = i. The electron density operator Pill(x) = 
&(x)*(x) can be bosonized as (l5||58[|59,61,62 


, JV dMx) « 
PlLL( ^= 


2nNx 


- 2ip(x) - 2f(x) 


(15) 
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with f(x) defined as in Eq. (11) and 

tp(x) sin (~J~) ~ bn ) ■ 

n>0 y ' U 


( i ) the density-density correlation function; 

(ii) the pair correlation function; 

(in) the probability density of finding two electrons at 
distance x. 


The zero-temperature average electron density P\ll(x) = (*) The density-density correlation function is 

(N\pill(x)\N) (with | N) the ground state of the ILL with 

N electrons) is d 1LL (N,x,y) = (N\p 1LL (x)p 1LL (y)\N). 


(19) 


Pill(x) 

K(x) 


N N 

l ~ l cos 


2Nnx 

L 


-2 /(*) 


K(x) (17) 


sinh(^) 


1 9 


\J sinh 2 (^) +sin 2 (?f) 


(18) 


and is shown in Fig. [2] for different interactions. 

Even in the noninteracting case g = 1, Fig. [2](a), the 




Fig. 2: Plot of Pill{x) (units 1/L) as a function of x (units L) 
for IV = 10 and different interaction strengths: (a) g = 1, (b) 
g = 0.5, (c) g = 0.1. In both panels a = L/(10n), T = 0. 


electron density exhibits N peaks. Increasing interactions 
only results in a continuous enhancement of the peak- 
to-valley ratio, Figs. [2jb,c) . Indeed, both the finite-size 
Friedel oscillations and Wigner oscillations have the same 
wavelength ~ L/N. Studying the density alone, it is 
therefore impossible to assess if a critical value g c (N) 
exists such that for g < <? c (iV) interactions overcome finite 
size effects. 

To detect the crossover towards a strongly interacting 
regime and finally to the formation of a strongly corre¬ 
lated Wigner molecule, which is the task of this Letter, 
we now compare three different correlation functions, 
namely 


Although its calculation can be carried out analytically, 
the explicit result is omitted for simplicity. A plot of 



Fig. 3: Plot of <1ill(N, x, y) (units 1/L 2 ) as a function of x 
(units L) for different interaction strength (a) g = 1; (b) g = 
0.2. In both panels N — 10, y = 0.65L, and a = L/(10n), 
T = 0. 


diLL(N,x,y) is shown in Fig. [3j We chose for y the 
location of one of the maxima in the electron density, in 
order to maximize the contrast. The function exhibits 
N distinct peaks, the auto-correlation one at x = y 
being dominant, especially for the non-interacting case 
shown in Fig. [3](a) . In the strongly interacting regime the 
only noticeable effect is an increase of the peak-to-valley 
ratio as seen in Fig. §b). Thus, this quantity behaves 
qualitatively in the same way as the electron density 
and no clear onset of an interaction-induced crossover 
towards the correlated Wigner regime for can be detected. 


(ii) The pair correlation function giLL(N,x,y) is [57 


giLL(N,x,y) 


(N\&(x)&(y)V(y)V(x)\N) 

Pill(x)pill(v) 


( 20 ) 


and has often been employed for characterizing strongly 
correlated systems [57. A plot of giLL(N,x,y) for 
N = 10 and y located at a maximum of the electron 
density is shown in Fig. [3] The most notable feature is 
the presence of a Pauli hole, namely the sharp collapse for 
\x — y\ <C L/N due to the Pauli exclusion principle. For 
enhancing the visibility of the oscillations we set the scale 
to a restricted region, the Pauli hole being unimportant 
for our aims. Both the noninteracting case (Fig. [4[a)) 
and the strongly interacting (Fig. |4](b)) have the same 
number of peaks, namely N — 1. Again, the effect of 
interactions is to enhance the peak-to-valley ratio with 
no additional feature emerging. We can therefore rule 
out also this quantity as a clear-cut estimator of the 
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Fig. 4: Plot of giLL(N,x,y) as a function of x (units L) for 
different interaction strengths (a) g = 1; (b) g = 0.2. In both 
panels N = 10, y = 0.65 L, and a = L/(10ty), T = 0. 


emergence of Wigner correlations in a ILL. 

(in) The probability density Pill 0*0 of finding two elec¬ 
trons at distance x has been recently proposed as a tool 
to detect Wigner correlations in ID [27|. It is defined as 


Pill (IV, x) = 


JZo d v( N \H X ’V)\ N ) 

N(N-l) 


( 21 ) 


where 


h(x, y) = ’P t (y + |) ¥ (y - |) ^ (y - |) ^ (y + |) . 

( 22 ) 

While the quantum average can be analytically performed, 
the integration is carried out numerically. Results for 


H 

-5 

cC 


-5 



of Pill(N,x) as x is increased is due to the reduction of 
phase space in the integration. The absence of correlations 
in Pill(N,x) confirms the liquid-like structure of the dot 
state for the noninteracting case. Increasing interactions 
to mild values, for g = 0.8 shown in Fig. [5](b), signatures 
of correlations begin to develop with Pill (IV, x) showing 
three distinct maxima for x < 1/2. Notice that due to 
the Pauli hole, the maximum number of maxima expected 
for Pill(N,x) is N — 1. This suggests to define g c (N) 
such that, for g < g c (N) one has N — 1 distinct maxima 
in P\ll(N,x). We have performed extensive numerical 
scans of Pill (IV, x) for 5 < N < 30 and found that the 
crossover is almost insensitive to the number of particles 
0.45 < <7 C (IV ) < 0.5. Indeed, Figs. [5jc,d) show Pill(N,x) 
for g < <7 C (10) ~ 0.46. Clearly, 9 distinct peaks located at 
x = x n ~ nL/ 10 (1 < n < 9) are present, which confirm 
the physical picture of an ordered molecular state. 

The origin of the better performance of Pill (IV, x) with 
respect to giLL(N,x,y) and cIill(N, x,y) lies in the inte¬ 
gration over the coordinate of the reference position with 
respect to which the distance x between correlated elec¬ 
trons is measured. Indeed, both in <?ill(IV, x, y) and in 
diLL(N,x,y) an unavoidable background of the uncorre¬ 
lated density oscillations remains even in the weakly in¬ 
teracting regime. Such a background is washed by the 
integration and genuine electronic correlations emerge. 

Conclusions. — We have studied the formation of a 
correlated Wigner molecule in a one-dimensional finite- 
size system of N electrons, described by a one-channel 
Luttinger liquid. Both Friedel and Wigner oscillations of 
the electron density are characterized by N distinct peaks 
located at the same positions. Thus, the density does 
not bring information on the crossover from finite size 
to correlation effects as interactions increase. We have 
compared three different tools to investigate correlations 
among the electrons: the density-density correlation func¬ 
tion, the pair correlation function, and the density prob¬ 
ability of finding two electrons at a given distance. We 
have shown that the latter is the best detector of the 
onset of Wigner crystallization among the three consid¬ 
ered in this paper. In a large range of particle numbers 
5 < N < 30, the critical value of the Luttinger parameter 
q r (N) below which the correlations appears is in the range 
0.45 < g c (N) < 0.5. 


Fig. 5: Plot of Pill(N, x) (units 1 /L) as a function of x (units 
L) and different interaction strength (a) g = 1; (b) g = 0.8; (c) 
g = 0.45; (d) g = 0.25. In all panels N = 10 and a = L/( 107r), 
T = 0. 
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IV = 10 are shown in Fig.[5]for x > a (for distances shorter 
than the cutoff a the calculation is not reliable). For zero 
interactions (g = 1), Fig. [5]) a), Pill(IV, x) exhibits a Pauli 
hole (whose precise form cannot be evaluated with our ap¬ 
proach for x < a) but is otherwise almost featureless with 
only a hint oscillations for x <C 1. The overall decrease 
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